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Abstract 



The decay Kl — ► "fw is analyzed within the standard model. Short- 
distance contributions are found to dominate, yielding a branching ratio of ~ 
0.7 xlO" 11 . We examine the possibility that non-standard model effects might 
be observed at higher rates. As an example, we calculate the branching ratio 
for the process mediated by neutral horizontal gauge bosons. Notwithstanding 
some model uncertainties, experimental limits for lepton number violation 
constrains these contributions to ^ 2 x 10~ n . 

(Accepted for publication in Physical Review D) 



1 



Typeset using REVTeX 



I. INTRODUCTION 



Rare kaon decays continue to be an area of active research |TJ. Within the standard 
model, rare decays are governed by loop level processes and thus directly probe the quan- 
tum structure of the theory. Experimental searches serve several complementary purposes. 
Improved limits help to constrain standard-model parameters, while discovery of clear dis- 
crepancies with standard-model predictions would signal new physics. Proposed extensions 
to the standard model are similarly constrained. 

Of particular interest are rare decays which are dominated by short-distance contribu- 
tions and hence are free of theoretical uncertainties from long-distance effects. A well-known 
example is K + — > tt + vv This decay is expected to give a clean determination of 

the Kobayashi-Maskawa matrix element Vtd, provided the top quark mass is well-known. 
Another example is Kl —>■ is^vv [@||, which is of special interest because its dominant 
contribution is expected to arise from direct CP violation. 

A rare kaon decay which has not been well-studied is Kl —>■ ^vv. Although the reasons 
for its study are perhaps not as compelling as those for the previous two processes, it is just 
as important to test this decay for consistency with standard-model predictions; as noted 
above, rare processes with very low rates are fertile ground in which to look for non-standard 
effects. Furthermore, information about this process will be useful in background analysis 
of other rare-decay experiments. An early calculation by Ma and Okada || predicted a 
negligible branching ratio of ~ 10 -18 . Their calculation, however, predates the confirmation 
of the 3rd quark family and clearly needs to be updated. In Sec. |I| we examine some 
long-distance processes and find contributions up to 0.7 x 10~ 12 . Short-distance effects are 
considered in Sec. |TT|, and we find a branching ratio of 0.7 x 10~ n . Thus within the standard 
model K L — > ^vv also appears to be short- distance dominated. 

One of the defining characteristics of the standard model is the absence of tree-level 
flavor- changing neutral currents (FCNC). However, many proposed extensions (technicolor 
theories for example) predict such currents. Rare neutral current kaon decays provide an 
ideal window to look for such effects. In Sec. IV we consider possible FCNC contributions 
that could arise from neutral generation-changing gauge bosons. Our concluding remarks 
appear in Sec. [V| 



II. LONG-DISTANCE 

The two long-distance contributions we consider are shown in Fig. [I]. Momentum and 
polarization assignments are indicated on the graphs. Each contribution is treated separately 
below. Note that within the standard model the form of the Kl — > decay amplitude 
is tightly constrained by Lorentz and electromagnetic gauge invariance. Assuming CP 
conservation, the only allowed form is 

e^(u^P L v)Ke pPa , (1) 

where u and v are neutrino spinors, and Pl — (1 — 7 5 )/2. 
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A. Kl — ► P —> ")w 

We first consider the long-distance diagrams shown in Fig. (HD(a). They involve a weak 
transition of the kaon to a virtual pseudoscalar meson P = Tr°,T],T}'. The corresponding 
amplitude A can be written in the form 

, y/2hpv\H\7r )(-K°\H\K } ^ A , 

A — -L-±l U ' / i + ^ + ^O, 2 

assuming { _P | | i^ ) = (Pl'^l-S' ). The quantities A v and A^/ represent the contributions of 
the r] and rj' diagrams relative to the pion graph amplitude. Their kaon transition elements, 
for example (t]\T-C\K°), can be related to the corresponding tt° transition element using SU(3) 
flavor symmetry. To account for interference effects, we carefully treat symmetry breaking 
0. We use the mixing convention 

\rj) = cos6» \r] 8 ) - sin6> \rj ), (3a) 
\rf) = sin 8 \r] 8 ) + cos 9 \rj ), (3b) 

where r] s and rjo denote the usual SU(3) octet and singlet states. The kaon transitions then 
satisfy, 



(n \H\K } ^3 

(vo\n\K°) 

(n°\n\K°) 




(4b) 



where ( and p parameterize SU(3) breaking. Current algebra gives (it \H\K°) ~ 1 x 10~ 7 m 2 K 
0. Chiral perturbation theory to one loop predicts ( ~ 0.17 @. The singlet parameter p 
and mixing angle are fit to experiment. We will take p = 0.83 and 9 = —20.6° ||. 

The pseudoscalar transitions (juv\7i\P) can be described by constituent quark triangle 
diagrams H], shown in Fig. The calculation mimics the well-known triangle graph 
approach to n° — > 77 0]. We obtain 



( TI /P \H\P) = ^-I p£ ^(u^P L v)k u e pPa , (5) 

where the pseudoscalar and quark loop dependencies have been collected into the coefficients 
I p. It is convenient to express Ip in terms of the mass ratios £ q = p 2 /m 2 for constituent 
quark flavor q. Then 

Ip = 2N c J2Q q 9 -^ C -^ farcsiny^) 2 , (6) 

provided < £ 9 < 4. Here iV c = 3 is the number of colors, Q g are the quark charges in units 
of e, and cy q are the standard-model weak vector coupling coefficients. The effective meson- 
quark couplings gp q q are shown in Table |. The parameters /§ = 1.25 f n and fo = 1.04f n , 
where f T = 132 MeV, account for SU(3) breaking in the rj 8 and rj amplitudes [Oj. For the 
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constituent quark masses, we choose the representative values m u = = 0.31 GeV and 
m s = 0.46 GeV. The resulting coefficient values are I w o = 0.25, I m = 0.51, and I m = 5.7 
GeV- 1 . 

In terms of Eqs. @, (f|) and (||), the quantities A v and A v i are 



^— 1=- cos 9 + \-p sin 9], (7a) 
V3^ /' 



V 2 2 \ r 1 t 

Inserting the numerical values for Ip gives = —6.9 and A v > = 10.6. Integrating the 
amplitude A over phase space gives the corresponding partial rate V for a single neutrino 
family 

1 " 60(2^ (m^-m^ (1 + A " + ^ ' (§) 

where a = e 2 /47r is the fine structure constant. For three neutrino families, we find using 
Eq. (Bp a branching ratio of 0.8 x 10 -17 . 



B. K L -> if* (7) z^P 



The second long-distance contribution we consider is shown in Fig. (P(b). The dia- 
gram involves a radiative transition of the kaon to an intermediate K*. The corresponding 
amplitude A takes the form 



.4 



({vv\H\K*°) + (vi?\H\K* )) ■ (K*°y\H\K ) 



m 



(9) 



K* j 



assuming (K* 0r )\H\K°) = (K*°y[H\K°). To model the photon vertex, we normalize the 
amplitude predicted by the vector meson dominance hypothesis to the experimental rate for 



K* -> K"y. The amplitude is JT2 



A K *^ K o 7 = ege lxupa q^E v k p e a , 
where g is the effective vertex coupling. The corresponding decay rate is 



a\g\ 2 m 3 K , 
24 



m 



K" 



m 



(10) 



(11) 



The measured branching ratio for K* — > K 0r y is 2.3 x 10 3 fL3 |. This implies \g\ = 1.3 
GeV" 1 . 



4 



The K* — > z/z/ transition is modelled as a short- distance process. In the standard model, 
the quark-level interaction ds — > vv proceeds at the one-loop level. As indicated in Fig. || 
both box and effective dsZ vertex diagrams contribute. The corresponding effective four- 
fermion interaction is [[| 

Hi = 2V2G FX Di s^P L d v^PLVi + h.c, (12) 

where \ = Ci/An sin 2 9\y and i = e,fi,r denotes the neutrino family. The quantity Di is a 
function of the appropriate lepton mass, the up-type quark masses, and several elements of 
the Kobayashi-Maskawa matrix. Our determinations for Di are described in the Appendix. 
The K* transitions determined by Hi contain the hadronic matrix element (0|s7 M P£,d|if*) = 
f K ,E p . We will take f K . = f p = (0.2)m 2 



Eqs. (0) and dl2|) determine all the transition elements appearing in Eq. The 
resulting amplitudes A{ for a single neutrino family are 

A = ^GF^ReA ^^^^^^^ (13) 
q — m K * 

The corresponding partial decay rates are 

_ a\g\ 2 {G F xfK*^Di) 2 m] { 

* 127T 2 

where x = m 2 K */m 2 K and F(x) is a kinematic function. Specifically. 



F(x), (14) 



17 x — 1 

F(x) = — - + 7x - Ax 2 - (x - l) 2 (4x - 1) In — — . 

Summing over the three neutrino families gives T = 8.5 x 10~ 30 GeV, or a branching ratio 
of 0.7 x 1(T 12 . 



III. SHORT-DISTANCE 

For Ki — > 7z/z/, the relevant short-distance interaction is ds — > ^vv. To classify the 
lowest order diagrams, it is convenient to start with the graphs for ds — > vv. These were 
shown earlier in Fig. |3|. Adding a photon to any charged line creates a graph for ds — > '-/vv. 
Provided the standard-model WWZj vertex is included, this method generates the full 
set of lowest order diagrams. The short-distance contributions thus naturally split into box 
(Fig. |) and effective dsZ type (Fig. ^|) diagrams. Each class further divides into one-particle 
irreducible and one-quark reducible processes. 

Ma and Okada || only considered the leading contribution to the effective dsZ type 
diagrams: the one-particle irreducible graphs containing a single W boson f2|. Assuming 
just two quark families, they found the following short-distance amplitude for Kl — > ^vv: 

V2eC? F sin^cos 9^J K ^ ( ^p Lv)k ^ (15) 
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where 6 C is the Cabbibo angle. Their result was dominated by the graphs containing an 
intermediate u quark. The heavier charm quark contributions were shown to be small. The 
analysis, however, did not take into account the possibility of an intermediate quark mass 
larger than the W mass. 

We have extended the Ma and Okada calculation to the case of three fermion families, 
and a large top quark mass. In the large mass approximation, the above amplitude would 
be modified by the factor 



P-k V td V t * s 
Qrri^r sin 6 C cos 9 C 



21nz 1 



[z-lY (z-1? zz-1) 



where z = m^/m^y. For m t ~ 174 GeV fl~5 |, this factor actually reduces the corresponding 



rate by several orders of magnitude. Thus Ma and Okada's conclusion that the effective dsZ 
diagrams are dominated by the lightest quark is still valid. 

Ma and Okada, however, did not consider the box diagram contributions. Our expecta- 
tion is that the major short- distance contributions in fact come from the one-quark reducible 
box diagrams in Fig. |](b). 

In order to later treat non-standard model effects, we approach the calculation in a rela- 
tively general manner. As indicated in Fig. || we consider all one-quark reducible processes 
that can be viewed as a photon attached to some effective ds — > VjVu interaction. The form 
of the four-fermion Hamiltonian is taken to be 

H = (3 jk sD^d UjN^ k + h.c, (16) 

where = 7^(vd — aols) and = r y fl (yN — a^vTs) are arbitrary mixtures of vector 
and axial-vector components. The standard-model effective four-fermion Hamiltonian was 
discussed previously in Eq. (|P2p. In this case, the fermion currents are both left-handed, 
lepton number conservation implies (3j^k = 0, and (3jj = 2\^2GfxDj- Note that by using 
these coefficients, several effective dsZ type graphs will be included along with the one-quark 
reducible box diagrams. 

Using the interaction H, the quark-level amplitude generated by the processes in Fig. |6| 

is 



Pd ■ k 
Ps ■ k 



UdUjN^Vk, (17) 



where and m s are quark current masses. It remains to extract the corresponding kaon 
decay amplitude. For a reasonable estimate we project out the "pseudoscalar" content of A 



with the operation [16 



M = Cj2u s i 5 v d A, (18) 

where the sum is over the quark spins and C is a normalization constant. Only terms that 
conserve CP are retained. In addition a simple static constituent quark model is used to 
evaluate the quark momenta. The resulting projected amplitude for K° — > 7iA,z/ fc is 
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M = ^f{m s - m i )v D e>w"'(p j N l tV k )k p (- v p a (19) 

where and m s are now constituent quark masses. The K° amplitude is identical except 
that P* k appears. Notice that only the vector part of the quark current contributes. 

The normalization constant C is determined by applying the projection to the de- 
cay K + — ► l + V[. Comparison with the definition of the kaon decay constant gives 
C = fKm K /8m s m d V* s . 

After combining the K° and K° amplitudes, and integrating over phase space, we obtain 
the corresponding K L — > jisv decay rate 

' " (v^ + a^HTr[(/5 + /3t) 2 ]%^, (20) 



2(72tt) 2V 1V ' ' - 1 \ Vui 

where p = m^/m^ — m d /m s . Inserting the standard-model couplings and coefficient matrix 
P gives 

r^BHeA)^. ( 21 ) 

Our numerical evaluation of the Inami and Lim function Di is discussed in the Appendix. 
Estimates for the constituent quark masses vary substantially. We will take p = 0.8 ± 0.2 



to account for some of this uncertainty. The remaining parameters are all well-known [|13 
The resulting branching ratio is 0.71q;4 x 10~ n . 



IV. HORIZONTAL GAUGE BOSONS 

Since K L — ► ^vv is highly suppressed in the standard model, observation of this decay at 
a higher rate would signal new physics. Therefore, it is of interest to look at the predictions 
of extended models. One possibility is the contribution of flavor- changing neutral currents 
(FCNC). Here we consider the effects of neutral generation-changing gauge bosons and, in 



particular, we follow the general formalism adopted by Cahn and Harari [|17| . We have 
extended this formalism to three generations, as follows. 

The general group structure is taken to be G x H, where G acts within a generation and 
if is a "horizontal" gauge group containing neutral generation-changing bosons. The group 
G contains the standard electroweak algebra SU(2)^ x U(l). The group H is taken to be 
SU(2)#. We will further assume that the three horizontal gauge bosons are degenerate in 
mass. 

The basis of particle states is constructed as follows. Let L°, L\ and L° denote the "prim- 
itive" electron, muon, and tau. For the present, it is convenient to ignore the distinction 
between left and right components. These primitive states are eigenstates of both $\]{2) w 
and SU(2)# with eigenvalues T 3 W = —1/2 and T 3 H = —1,0, 1 respectively. Similarly, let the 
triplets Uf and denote the primitive states of the neutrinos, up- type, and down-type 
quarks. 

The set of mass eigenstates is assumed to be related to the primitive states by a unitary 
transformation. Specifically, 
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' N ' 




\U N 




' N° ' 




L 




U L 




L° 


F = 


U 




u u 




U° 




D _ 




U D 




D° 



= UF° 



In the primitive basis, the horizontal currents have the simple representation p°T H F° where 
T H = IL4 ® t and the T\ are the generators for the 3- dimensional representation of SU(2). 
For convenience the Lorentz structure of the currents has been temporarily suppressed. 
The horizontal gauge bosons thus presumably give rise to the following low-energy effective 
interaction 



n 



2m 2 H 



(FUrU ] Ff 



(22) 



where gn is the gauge coupling constant and tuh is the mass of the gauge bosons. 

For the decay Kl — > 7^, the relevant interactions in Ti are those that connect the D 
and N triplets. These can be written as 



U DN = ^MrDrD ■ NU n U d ^tU d U n ^N, 



ni 



(23) 



H 



by making the replacement t — > U d ^tU d . Furthermore, keeping just the terms involving 
the d to s quark transitions gives 



n 



ds 



V2m 2 H 



(U N U m TM D U m ) jk sd Uju k + h.c. 



(24) 



Comparison of this expression with the general Hamiltonian given in Eq. (jig) determines 
the coefficients /3j k that appeared in our short- distance analysis. Thus the decay rate given 
in Eq. (p0|) can be used directly. Note that the unitary mixing matrices do not contribute 
to the trace. The resulting Kl —>■ ^vv rate is 



r 



2a 



(72tt) ; 



(m 2 j_ ^>,2 9h fKP 2m K 
\ N d N! y D~ 



m 



H 



V, 



(25) 



The neutral generation-changing bosons that we are considering also mediate the lepton- 
number violating processes Kl — ► e/i and K + —>■ 7c + efi. The stringent experimental limits on 
these processes help to constrain the horizontal gauge boson couplings and mass. The process 
Kl —>■ efi, however, is only sensitive to the axial-vector component of the quark currents. 
Thus the semi-leptonic decay K + — > TT + e[i, which is sensitive to the vector component, is 
more convenient as a comparison process. 

The relevant quark-level process for K + — > 7i + e~fi + is ds — * e~fi + . Repeating the 
steps for the DL interactions that were used for the DN interactions, one finds that the 
appropriate term in the Hamiltonian is 



H 



sde/i 



(26) 
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where V = U L U D ^ governs the mixing between lepton and down-type quark states. Normal- 
izing the corresponding decay rate to the standard-model decay K + — > 7r°i/ /1/ u + we find the 
ratio 

F „4 h; v* _i_ V V* I 2 

D _ 1 ?r+e-/i+ / 2 , 2\ 2 SH l^e^^MT" 1 " "ee | 

K = t; = l v L + a LJ v D" 



v ^ Ly "mf, G£|V U . 



The present experimental limit is i? < 6.6 x 10 [13 



The Kl — > 7z/z/ rate appearing in Eq. fl25|) can be written in terms of -R as 



In order to obtain a reasonable estimate, we ignore the remaining dependence on horizontal 
coupling coefficients and mixing angles. We then obtain for Kl — > ^vv a branching ratio 
limit of < 2 x 10 -11 . 



V. CONCLUDING REMARKS 

Within the standard model, the decay Kl — > ~jvv is indeed very rare, but not nearly to 
the extent predicted by earlier calculations 0. We have found long-distance contributions 
with branching ratios up to 0.7 x 10~ 12 and short-distance contributions of 0.7i !4 x 10~ n . 
We also considered the possible effects of neutral generation-changing gauge bosons. Present 
experimental limits for lepton violation suggest that these contributions are limited to < 
2 x 10 -11 . 
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APPENDIX: 

This appendix discusses the numerical evaluation of the Inami-Lim function Di || . Val- 
ues for quantities not specified in the text are taken from the 1994 Review of Particle 
Properties |T3 . 



Constraints on the KM Matrix 

Since Di depends heavily upon the elements of the Kobayashi-Maskawa (KM) matrix, 
we first update their values ||18| . In the Wolfenstein convention |I9]| , to order A 3 , the KM 
matrix V takes the form 
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V 



1 - A 2 /2 A A\ 3 (p-iri) 

-A 1 - A 2 /2 AX 2 

AX\l- p-iri) -AX 2 1 



(Al) 



The measured element magnitudes give A = 0.221 ±0.002 and A = 0.82±0.10. Furthermore, 
the ratio \V u b/V c b\ describes a circle in the (p,rj) plane, 



2 i 2 
p +7] 



(0.36±0.09) 2 . 



(A2) 



Two additional constraints on p and r], as reviewed below, are obtained from measurements 
of K°-K° and B° d -B° d mixing. 

In the neutral meson system, mixing is characterized by the parameter x^. The 
dominant processes are box diagrams containing top quarks. One obtains p0 | 



G 2 F m 2 
6tt 2 



W r2 



f B B B m B TB\E(x t )\r) B \V tb V td \ 



(A3) 



where the function E(x) denotes E(x,x), and || 



E(x,y) 



xy_ 
4 



V 2 - 8y + 4 
(y-l) 2 (x-y) 
3 



In 2/ 



2(x-l)(y-l) 



v) 



(A4) 



The arguments x q = rd\ |/m^ where m 9 are the quark masses. In the Wolfenstein convention, 
Eq. ( |A3| ) describes a second circle in the (p, 77) plane, 



V 



;i.l9±0.29) ; 



(A5) 



where we used f B = 0.140±0.025 GeV and B B = 0.85±0.10 |^T| for the hadronic parameters, 
m t = 174 ± 16 GeV [HJ for the top mass, and r] B = 0.85 ±0.05 [51] for the QCD correction. 

In the neutral kaon system, mixing is characterized by the CP violation parameter e. 
Both charm and top quark box diagrams contribute. One finds [2^] 

1 1 G 2 F m 2 w f 2 K B K m K 



I2n 2 ^Am K 



E(x i ,x j )T]ijImA i A 



i=c,t 
j=c,t 



where A q denotes the product V q dV* s . Note that when evaluating the A q , the KM matrix 
elements must be expanded to order A 5 . Using B K = 0.7 ±0.2 and fx = 0.160 GeV [^3|] for 
the hadronic factors gives the hyperbola 

(7.3 ±2.1) x 10~ 4 



(A6) 
0.76, 



The one-sigma contours for the constraints (|A~2|), (|A5|) , and ( |A6| ) are plotted in Fig. |^. 
The product of these distributions determines our expected region for p and 77. Contours 
for this region are shown in Fig. Rl 



V A 2 \ Vcc E(x c ) - Vct E(x c ,x t ) - A 2 A 4 (1 - p)vttE(x t )\' 

For calculations, the charm mass was taken to be m c = 1.5±0.1 GeV, and we used r\ c 
rj ct = 0.36, and r\ tt = 0.61 [123|] for the QCD corrections. 
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Evaluation of Di 

In the Wolfenstein convention for the KM matrix, the function D,, j 6 {e, fi, r}, becomes 
A ~ -Xr] c D(x c , yi) - X 5 r] t A 2 (l - p - ir])D(x t , y^, 
where 7] c and rj t are QCD corrections and 



D(x,y) 



x 



y v 



+2-3 



x-y\y-l 
y + 2 



Iny + (x <-> y) 
x 2 - 2x + 4 



(A7) 



+ 



lnx 



For the standard- model masses j/< 1, and y x except when y = y T and x = x c . It follows 
that -Dj is relatively insensitive to i. We again take m c = 1.5 ± 0.1 GeV and m t = 174 ± 16 
GeV, and set r] c = 0.7 and r\ t — 1 [24]. Then using the Wolfenstein parameter values 
determined above, we find 



ReD 



e,/i 

ReD T 



-1.8+°;^ x 10~ 3 , 

0.6 ^ 1 n -3 

0. 

2.0 ^ i n- I 



-1.6^ x 10 



(A8) 



ImA = 3.01^3 x 10" 

The listed values are the modes and the uncertainties enclose 68.3% of the distributions. 
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FIGURES 



FIG. 1. Long-distance contributions. 

FIG. 2. Triangle diagram for the pseudoscalar decays P — > ^vv. The effect of the diagram 
with the photon and Z lines interchanged is to double the weak vector couplings cy q and cancel 
the weak axial- vector couplings CA q - 

FIG. 3. Diagrams for ds — > z^za The four effective dsZ vertex diagrams can be classified as 
either (a) one-particle irreducible or (b) one-quark reducible. Only one of each pair is drawn, (c) 
The box diagram. 

FIG. 4. The box diagrams for ds — > ^vv can be generated from the box diagram for ds — > z^P 
by photon emission from all charged lines. Shown here are examples of (a) one-particle irreducible 
and (b) one-quark reducible graphs. 

FIG. 5. Effective dsZ diagrams for ds — ► 7^z/ can be generated from the effective dsZ diagrams 
for ds — > z^z/ by photon emission from all charged lines. Shown here are examples of (a) one-particle 
irreducible and (b) one-quark reducible graphs. 

FIG. 6. Effective short-distance processes. 

FIG. 7. One-sigma constraint curves for p and 77. The dashed and dotted lines respectively are 
the V u b/V c b and B mixing circles. The solid lines are generated from the K mixing hyperbola by 
a Monte Carlo treatment of the uncertainties. 

FIG. 8. Expected (p, rj) region. Dashed contour encloses 68.3% of the distribution; solid contour 
encloses 90%. 
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TABLES 



TABLE I. Effective meson-quark couplings. 



p 


9Puu/m u 


9Pdd/ m d 


gpss/m s 


Vo 


(2) 1/2 /^ 

(2/3) 1 / 2 /^ 1 
2(3)~ 1 /2 /(r i 


(2/3) 1 / 2 / 8 - 1 
2(3)-V2 /o - 1 




-(8/3)V2/ 8 -i 
2(3)- 1 / 2 / " 1 
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(c) 
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